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1. Introduction 



In this paper, we obtain several expansions for ({s) involving a se- 
quence of polynomials in s, denoted by ak{s). These polynomials can 
be regarded as a generalization of Stirling numbers of the first kind and 
our identities extend some series expansions for the zeta function that 
are known for integer values of s. The expansions also give a different 
approach to the analytic continuation of the Riemann zeta function. 

The inspiration for our formulas comes from Renter's short note 
in the May 1999 Monthly [K] where he derives a formula for Euler's 
constant 7 which can be regarded as the s ^ case of (11.81) . after 
subtracting 1/s from both sides. 

We start with Riemann's formula 



roo s-1 

Jo e — i 



dx 



and substitute t 



Let 



(if 

/ (-log(l-t))-i^, ^s>l. (1.2) 
Jo ^ 
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and consider the Taylor expansion 

oo 


We will derive the following recursion: 

ao{s) = 1 , 

s-1 

"i(^) = > 

_ 1 ' a,{s)j(k + e + s{2k + 2-j)) 

"'=+1^^^ - ;t(A; + l)(A; + 2)^ ^k-j + l){k-j + 2) ' 

(1.5) 

and prove the following theorem concerning ak{s): 

Theorem 1.1. For k > 1, ak{s)/{s — 1) is a polynomial in s with 
positive rational coefficients, and ak{s) satisfies: 

(l + log(A; + l))N+^ 
< Cs , (1.6) 



where 

Is- II 

\s\ + 

We then have the following identities: 
Theorem 1.2. 



^(|s| +2)21^1+1 . (1.7) 



rw^f^. (1.8) 

For positive integer X: 

r{s)c{s-x) = f:Ms)t(-^r \fk-f-i ^'-''^ 

k=0 j=l ■' 

t'o h '"is + k-j-2)...is + k-X-l) 

(1.11) 
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where S and E respectively denote the Stirling numbers of the second 
kind and the Eulerian numbers. We also have: 

oo oo oo / \ 

fc=0 fc=0 n=0 y ^ ^ ) 

with ^i{s + k) the trigamma function. 

Equation (11.91) is known in the case that s is a positive integer where 
the coefficients ak{s) can be expressed in terms of the Stirhng numbers 
of the ffist kind. Jordan [Jl Sec. 68, (11), pg 194] credits it to Stirhng. 
See also |RSj and [S]. This is described in Section [31 



2. Recursions and bound for ak{s) 
To study the q;a;(s)'s in greater detail consider 

jMt)'-') = {s-l)g{tr'Y.TT^''^ (2-1) 





so 



j^{9{tr')9{t) = {s- l)g{tr' E ^ • (2-2) 

On the other hand, differentiating (11.41) term by term, the left hand 
side above equals 

. oo X ^ oo ^ X 

E(fc + iK+i(^)^' Er-T • (2-3) 

\ ^ ^ ^ 

Equating coefficients of t^ in (12. 2p and (12. 3p . we get 

y ^'7'^"^-^^^'^ ={s-l)y aAs) ^~'^^ . (2.4) 
Using ao{s) = 1, we write this as 
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The first few aj(s)'s are listed below 



"2(5) 
0:3(3) 
04(5) 



1 

s - 1 



(2.6) 



— s H s H 

48 16 24 



s'' + 

384 384 



23 



19 \ 



s + 
576 720 J 



We also observe, from fl2.5p . that we can write 



k + 2 k + 



j-(g-l)(fc-j + l) 
A;-j + 2 



(2.7) 



Furthermore, substituting k — 1 for A; and moving the l.h.s. to the r.h.s., 
we have 

.-1 (2.8) 







k + l k 



A; - J + 1 



Subtracting {k + l)/{k + 2) times ([23]) from we get 



a,(s)j(A; + A;2 + s(2A; + 2-j)) 



A;(A; + l)(A; + 2) ^ (k - j + l){k - j + 2) 



k>l. 
(2.9) 



This last form for the recursion governing the a^'s is useful in that we 
easily obtain, inductively, the first part of Theorem 11.11 namely that 
ak{s)/{s — 1), for > 1, is a polynomial in s with positive rational 
coefficients. 

Next, we determine a bound for |Q;fc(s)|. Let Mg be the positive 
integer satisfying |s| + 1 < Mg < |s| + 2. From the positivity of the 
coefficients of ak{s)/{s — 1) we have 

ak{Ms + l) 



s-1 



< 



A; = 1,2, 



(2.10) 



Here we are using the fact that ak{s) /(s — 1) is a polynomial in s with 
positive coefficients and also Mg + 1 > \s\. Rewriting fl2.10p , and using 
|s| + 1 < Mg, we find 



\akis)\ < 



Mg 



-ak{Mg + 1) < 



|s| + 1 



ak{Mg + 1) . 



(2.11) 
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Now, 

Lemma 2.1. 

\MM + 1) I < ' M = 1, 2, 3, . 



Proof. 

ak{M + 1) is the coefficient of t'' in g{t)'^'^ = ( 1 + ^ + V + 



t ^^ 



2 3 

When M = 1 the lemma is easily verified. Now assume that the lemma 
has been verified for M and consider the M+1 case. Writing g{t)'^~^^ = 
g{tY'g{t) we get 

k 

ak{M + 2) = V ^-a,_^(M + 1) . (2.12) 
By our inductive hypothesis, 

MM + 2)1 < V _^M2"-'(l + l°s(^-n.+ l))"-' 

^ ^' - ^ m + 1 k-m+1 

M2^i-i 1 / 1 1 



^(l + log(A;-m + l))^^-i C 



k+2 ^ \m+l k — m + I 

m=0 

<^2-'(l + log(. + l))"-t;^. (2.13) 

m=0 

But 2 J2t=o + 1) < 2(1 + log(fc + 1)), and M/{k + 2) < (M + 
l)/{k + 1), thus the lemma is proven. □ 

Combining (12.111) with Lemma 12.11 and the assumption that Mg < 
|s| + 2 yields the bound in Theorem 11.11 

< ^ (|.| + 2)2M^Hl + log(t- + l))'-'^' ^ p,14) 

1 5 1 1 k -\- \ 

We conclude, by using the above bound to show that the sum below 
converges, that we may substitute (11.41) into (II. 2p and integrate the 
power series term by term to obtain 

^ s + fc - 1 

Away from its poles, the above is, in bounded sets, a uniformly con- 
vergent sum of analytic functions. Thus, while we started in (II. ip with 
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> 1 , formula (12.151) gives the meromorphic continuation of the left 
hand side. 

Similarly, applying the same change of variable to the integral defin- 
ing the Gamma function, 



oo 



e dx 

(2.16) 



f {-\og{l ~t)y-^dt , > 

Jo 



gives 



oo 



rw = X:^. (2,17) 

3. Connection to Stirling numbers of the first kind 

The coefficients ak{s) defined by (11.41) are related to Stirling numbers 
of the first kind through the series expansion, for integer m > 0, 



oo 



log(l + tr = m! ^ ^^^t"- (3.1) 



Thus, when s is a positive integer. 



a,{s) = i-^) \Xk-[)& + k-l,s-l), (3.2) 
and equation (I2.15P becomes 

As mentioned in the introduction, the latter formula, valid for positive 
integer s, is known and (12.151) may be regarded as a generalization of 
(13.31) to s G C. The expression (13. 3p can be combined with identities 
that relate the Stirling numbers to the harmonic numbers, and from 
this various identities for zeta evaluated at positive integer values may 
be deduced. See for example Section 4 of [RS] or [S]. 

4. Connection to Stirling numbers of the second kind 

Other related expansions of ({s) are possible. For example, instead 
of (11.11) . consider 

r(s)C(s-A)=/ x'-^VnV"^rfx, 3?s-A>l. (4.1) 
Jo 1 
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Substituting t = 1 — e^^, we can rewrite the above as 

When A is a positive integer, we can evaluate the sum over n by 
repeatedlv multiplying by (1 — t) and applying —d/dt to the geometric 
series Ef (1 " ^)""^ = V^- 

oo 



Therefore, 



1 1/f 

4 24/^5-36/^^ + 14/^3-1/^2 
akis) 



r(.)C(.- 2) = ^ 

^ 

oo ^ 

r(.)C(.-3) = 

^ 

oo X 

n V 



s+k-3 s+k-2 



6 6 1 

+ 



s+k-4 s+k-3 s+k-2 



24 36 14 

+ 



s + A;-5 s + A;-4 s + k-3 s + k-2 



Next, we derive the general form for the above expressions and give a 
connection to Stirling numbers of the second kind: For positive integer 
A, 

00 A 

^n^(l - tr-' = ^(-l)'+^i!^(A, j)/t^'+' . (4.4) 

1 j=i 

We can verify this formula inductively, multiplying by (1 — t), applying 
—d/dt, and using the the recurrence relation for Stirling numbers of 
the second kind: 

S(X + 1, j) = S{\ j-l)+ jS(X, j) . (4.5) 
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Substituting into (14.21) and integrating gives the formula 

r(.)c(. - A) = f; M.) h-^}''' Jf^'t 1 («) 

fc=0 j=l 

of Theorem 11.21 

We can also go in the opposite direction, taking, for example, A = — 1 
in (g^D- Now, 

f:n-'(l-tr- = i2li. (4.7) 

1 

Using 

\'+''-^^^dt = ^^{s + k) , (4.8) 



where is the trigamma function, we get 



oo oo 



k=0 k=0 n=0 ^ ^ 

5. Connection to Eulerian numbers 

We can also expand the sum in (14.21) in terms of Eulerian numbers. 
For positive integer A, 

00 A— 1 

Y,n\i - tr-' = t-'-' E ^(^'^^(i - tf~^'' , (5.1) 

1 0=0 

where E{X,j) are the Eulerian numbers, satisfying the recursion: 

EiX + 1, j) = (j + l)E{\,j) + (A + 1 - j)^(A, J - 1) . (5.2) 

Formula (15.11) can, again, be verified inductively, multiplying by (1 — 
t), applying —d/dt, and manipulating slightly before using the above 
recursion. 
However, 

f ^^-^-(1 - tr^-'dt = ns+k-x-mx-j) 

T{s + k-j- 1) 
i.e. the Beta function, gives, after applying T{z + 1) = zT{z), 

(5.4) 
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6. Obtaining C(1-A), A = 1,2,3,.... 

We can apply formula f ll.lOp to obtain expressions for ^(1 — A), with 
A a positive integer. We first illustrate the technique for A = 1. 

By Theorem II. ![ ak{s), for A; > 1, is a polynomial in s divisible by 
s — 1. Thus, ak{l) = 0, k > 1. We also have ao(l) = 1. Therefore, 
substituting s = 1 into the first equation in (14. 3p . all but the first two 
terms vanish. The k = term gives —1, while the denominator of the 
k = 1 term cancels the s — 1 factor of ai{s) = (s — l)/2, and the k = 1 
term equals 1/2. We thus get 

C(0) = -1 + 1/2 = -1/2 . (6.1) 

The general situation is handled via equation (ll.lOp . As in the A = 1 
case, we get two kinds of contributions- from the k = term, and from 
the terms k = j, with j = 1, . . . , A. The latter terms are the ones for 
which the denominator s + k — j — 1 cancels the factor s — 1 of ak{s). 

The k = term produces, on simplifying, 

A 

i-irY.(-^y''(3-iy.s{x,j) . (6.2) 

But, the Stirling numbers of the second kind are defined by the expan- 
sion 

A 

x^ = Y,S{X,j)x{x-l)...{x-j + l) , A>1. (6.3) 

Thus, dividing by x, and setting x = shows that (16. 2p equals if 
A > 1, and —1 if A = 1. Therefore, the k = term only contributes 
when A = 1, and we denote its contribution by —6x, equal to or — 1 
according to whether A > 1 or equals 1. 

The terms k = j, with j = 1,2,..., A, contribute to (ll.lOp . on 
canceling the zero and pole at s = 1, the sum: 

A 

(-l)^^a,(l)'(-l)'=A:!^(A,A;) . (6.4) 

k=l 

Therefore , putting both contributions together: 

A 

C(l - A) = -6x + (-1)" Yl M^n-^)'k\S{\, k) . (6.5) 

k=l 

We thus need a formula for afc(l)'. We can differentiate (12.70 . and 
use the fact from Theorem 11.11 that (s — l)ak{s) has a double order 



10 



MICHAEL O. RUBINSTEIN 



zero at s = 1 when A; > 1, to get 

k 

This recursion allows us, in conjunction with fl6.5p . to evaluate (^(1 — A) 
for any positive integer A. 

Equation (16.51) can also be used to derive Euler's famous formula 
involving the Bernoulli numbers. To see the connection to Bernoulli 
numbers, we calculate the first few values of afc(l)' from the above 
recursion, starting with ai{l)' = 1/2, and list them in Tabled! 



k 










1 


1/2 


2 


5/24 


3 


1/8 


4 


251/2880 


5 


19/288 


6 


19087/362880 


7 


751/17280 


8 


1070017/29030400 


9 


2857/89600 


10 


26842253/958003200 


11 


434293/17418240 


12 


703604254357/31384184832000 


13 


8181904909/402361344000 


14 


1166309819657/62768369664000 


15 


5044289/295206912 


16 


8092989203533249/512189896458240000 


17 


5026792806787/342372925440000 


18 


12600467236042756559/919636959090769920000 


19 


69028763155644023/5377993912811520000 


20 


8136836498467582599787/674400436666564608000000 


Table 1. Values of 0:^(1)' 



Googling some of the larger numerators in this table, for example 
703604254357, immediately returns entries A002208 and A002657 from 
Sloane's Online Encyclopedia of Integer Sequences [Si]. These entries 
deal with the numerators of the Norlund numbers, and the following 
formula is stated: 'Numerator of integral of + l)...(x + n — 1) from 
to 1.' Comparing a few of our numbers to the values of these integrals, 
we find that the denominators in our case are off by a factor oi k k\ 
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from the denominators of this formula. We are thus led to surmise 
that: 

Lemma 6.1. The following formula holds for k > 0: 



1 



afc(l)' = J {x)kdx , (6.7) 
where 

{x)k = x{x + 1) . . .{x + k -1) . (6.8) 

Proof. We show that the r.h.s. in the lemma satisfies the same recursion 
as ak{iy. Replacing, in (16.61) . the a's by the r.h.s. of the lemma, 
rearranging and simplifying slightly, we wish to show that 

fc+i . . .1 , . 

x)jdx = -j-—^ , A; > 0. (6.9) 



^-.k-j + 2j\Jo " k + 2 

The l.h.s. above is the coefficient of z^'^'^ in the product: 



,=1 ^ -J^.r.Jo 



The first sum is the power series for — log(l — z), while the second sum 
equals the integral from x = to 1 of 



OO / 

[x 



(6-11) 



7 

which is the power series for (1 — z)^^' — 1. Integrating from to 1, we 
find that the product in (I6.10p equals 

- - ^' ( (.-i)iog(i-.) - = + - 

whose coefficient of is 1 - 1/(A; + 2) = (A; + l)/{k + 2). □ 

Substituting Lemma 16.11 into (16. 5p and rearranging summation and 
integration gives 

,1 A 



c(i - A) = -5. + (-1)^ / 

-^0 fc=l 



S{\k) 



{x)kdx . (6.13) 



k 

Moving the (— l)'^ into the {x)k and changing variables, u = —x, yields 

-6x + (-1)^ r V ^^%^m(u -l)...{u-k + l)du . (6.14) 
-^-1 k=i ^ 
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Next, apply the recursion (14. 5p to split the above sum over k into two 
sums. The second of these sums equals 

A-l 

^ 5(A - 1, k)u{u - 1) . . . (m - + 1) , (6.15) 

k=l 

which, by (16. 3p is u^~^ , if A > 1. If A = 1 it equals 0. Integrating from 
—1 to 0, it contributes to fl6.14p : 

-1/A ifA>l, 
^0, ifA = l. 

which we can write as 

-1/A + 5a. (6.16) 
Therefore, (16.140 has been simplified to 

((1-A) = -i/A+(-l)^ r ^ ^^^~^'^^~^\ {u-l) . . . {u-k+l)du . 

•' k=l 

(6.17) 

However, the sum above can be related to sums of powers, and hence 
to Bernoulli polynomials. Assume for now that m is a positive integer. 
First, we use the operator A/(m) = /(m + 1) — /(m), applied to 
m(m — 1) . . . (m — + 1) to get rid of the numerator above: 

(m+1) . . . (m— A;+2)— m(m— 1) . . . {m—k+l) = A;m(m— 1) . . . {m—k+2). 

(6.18) 

Dividing by k, summing over m = to m — 1, and telescoping yields 
the well known identity 



u(u — 1) . . . (u — k + 1) ^-4 , N / , X 

^ i — -A L = ^m{m-l)...{m-k + 2), k>l. 



m=0 



(6.19) 

(The m = term is needed ii k = 1.) Therefore, the sum in (16.170 
equals 

A u-1 

^ 5(A - 1, - 1) ^ m(m - 1) ... (m - A; + 2). (6.20) 

k=l m=0 

Rearranging the two sums, and using 

A 

^ 5(A - 1, A; - l)m(m - 1) . . . (m - A; + 2) = m^"^ , (6.21) 
k=i 
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the integrand in f l6.17p becomes, for positive integer u, 



u-l 



^m^^^ (6.22) 



m=0 

But, sums of powers can be expressed in terms of the BernouUi poly- 
nomials and numbers: 

X:m^-^ = ^?^M^, A>1, (6.23) 

giving a formula for the integrand which is valid for all u, and not just 
positive integer u, since both the integrand in (16.171) and the above are 
polynomials in u agreeing on infinitely many values. We thus have 

C(l-A) = -l/A + (-l)^y"° , A>1. (6.24) 

Furthermore, 

J Bx{u)du = (-1)^ (6.25) 

which can be obtained by substituting t = u + 1, using the difference 
equation for the Bernoulli polynomials, 

Mi±ii^^ = „^ A>1, (6.26) 
A 

and applying the third of the defining properties of the Bernoulli poly- 
nomials: 

Bo{t) = 1 

B',{t) = kBk-iit), k>l 

1 

Bkit)dt = 0, k>l. (6.27) 
Plugging (16.231) (16.251) into (16.241) therefore gives Euler's formula: 

C{l-X) = (-l)'+^^ , A>1. (6.28) 

7. Further properties of ak 

Because Bx = for odd A > 1, Euler's formula gives, as is well 
known, ({—m) = for positive even integer m. Hence, the left hand 
side of (I2.15P has poles at s = 1, 0, —1, —3, —5, —7, . . .. Therefore, to 
cancel the poles of the right hand side at s = —2, —4, —6, . . ., ak{s) 
must be divisible by s + A; — 1 when A; = 3, 5, 7, . . .. 
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Similarly, from the first equation in (14. 3p . ak{s) is divisible by s+k—2 
when = 1, 3, 5, 7, . . ., and from the second formula, 2ak+i{s) — ak{s) 
is divisible by s + A; — 2 when = 0, 2, 4, 6, . . .. 

By comparing the residue of (12.171) ats = — A;,A;>0, we get 

aki-k) = (-I)VA;!. (7.1) 

Next, by considering the residue of (12.151) at the poles s = —1, —3, —5, . . ., 
we have 

a2m+2{-2m-l)= , m>0, (7.2) 

[2m + 2)! 

and from first formula in (14. 3p . 

«2„^+2(-2m) = -(2m + 1) . , m > 0. (7.3) 

(2m + 2)1 
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